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Abstract 


The  situation  of  a  circular  wall  Jet  injecting 
into  a  channel  flowing  in  the  same  direction  is  inipo 
that  it  directly  parallels  some  of  the  methods  used 
effluent  disposal  in  marine  and  river  environments* 
Research  was  conducted  into  the  extent  and  the 
of  downstream  spread  of  the  jet  to  give  some  methods 
orediction  of  its  velocity  at  any  point  in  the  flow* 
Since  previous  research  has  shown  that  the  velo 
patterns  produced  by  jets  of  various  nozzle  shapes  w 
similar,  a  circular  jet  was  chosen  for  simplicity.  T 
downstream  velocity  profiles  were  found  to  be  simila 
been  predicted  for  free  jets  and  wall  jets  issuing  i 
stagnant  environment* 

Relationships  between  downstream  velocities  and 
excess  mo men  turn  thickness  ©  that  the  jet  posesses  we 
to  be  similar  for  the  various  velocity  ratios  invest 
and  the  theoretical  justification  for  the  use  of  thi 
is  shown  in  the  analytical  section* 

The  c orrespondence  between  the  velocity  at  any 
and  the  concentration  thet  results  from  an  injection 
pollutant  at  some  position  upstream  is  also  shown  so 
these  results  can  be  used  for  the  prediction  of  that 
concentration  at  any  downstream  section* 

Agreement  between  these  results  and  the  relatio 
their  theoretical  prediction  was  extremely  good* 
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I*  Introduction 


The  problem  of  effluent  disposal  in  marine  and  river 
environments  greatly  concerns  engineers  and  ecologists  in 
the  present  day.  With  the  rapid  expansion  of  cities,  towns 
and  small  communities  the  disposal  of  associated  domestic 
and  industrial  sewage  requires  considerable  thought  and 
planning  so  as  not  to  either  instigate  or  increase  harmful 
pollution  in  the  receiving  waters. 

Large  bodies  of  water  have  the  ability  to  absorb 
relatively  small  amounts  of  wastewater,  containing  suspended 
solids  and  biodegradable  products,  and  purify  it  with  no 
detrimental  effects.  It  is  when  relatively  large  quantities 
of  wastewater  are  introduced  to  a  water  body  that  severe 
local  and  global  pollution  may  result. 

It  is  therefore  clear  that  the  most  serious  effects  of 
this  pollution  will  be  experienced  in  inland  waterways  which 
must  carry  the  effluents  of  large  cities.  As  these  rivers 
and  lakes  often  provide  recreation,  and  also  water  supply 
for  the  city,  it  is  all  the  more  important  that  they  should 
be  kept  as  free  from  any  serious  pollution  as  is  reasonbly 
po  ssi b le  • 

The  main  criterion  when  judging  the  possible  effects  of 
wastewater  introduction  is  the  dilution  that  wastewater 
experiences  a  short  distance  from  the  point  of  injection.  It 
is  generally  accepted  that  a  dilution  of  one  part  of 
wastewater  to  one  hundred  parts  of  receiving  water  (  1:100) 
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The  disposal  of  wastewater 
presents  much  less  of  a  problem 
counterpart.  Vast  volumes  of 
the  wastewater  of  even  large  cities 
exist  enabling  the  outfall  to  be  placed 
water  surface  as  possible.  The  most 
would  be  if  the  effluent  never  reached  the  water 
however  with  the  density  difference  between  salt 
fresh  wastewater  the  effluent  possesses  buoyancy 
naturally  rise  towards  the  water  surface  unless 
density  gradients  exist. 

In  river  situations  the  density  difference 
pronounced  leading  to  effluent  plumes  and 
the  surface  of  the  water  much  more  slowly. 


Fig.  1-01  is  a  passive  plume  of  dyed  fresh 
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FIG.  1-01  SKETCH  OF  SINGLE  RISING  PLUME  OF  FRESH  WATER  IN 
A  STAGNANT  SALT  WATER  ENVIRONMENT. 
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stagnant  salt  water  environment.  The  fresh  water  is  passed 
through  a  horizontal  outlet  where  immediately  buoyancy 
forces  cause  it  to  rise  to  the  surface  and  spread  out  in  a 
layer.  Rawn  et  al  (1960)  analysed  the  dilution  that  occurred 
and  found  it  to  be  a  function  of  the  ratio  of  water  depth  to 
nozzle  diameter  and  Froude  number  at  the  outlet.  Below  a 
Reynolds  number  of  5000  this  parameter  also  becomes 
important;  however  almost  all  flows  will  be  well  above  this 
va lue • 

It  is  reasonable  to  assume  that  effluent  disposed  of 
into  rivers  will  also  rise  similarly  as  there  will  be  a 
certain  amount  of  buoyancy  due  to  the  excess  heat  of  the 
wastewater.  This  is  approximately  the  situation  that  occurs 
in  the  disposal  of  effluent  through  an  outfall  pipe  into  the 
ocean.  It  does  not*  however,  usually  produce  high  initial 
dilutions  and  is  only  used  if  the  water  is  deep,  and  away 
from  any  sites  of  recreational  activity  as  it  takes  longer 
for  the  wastewater  to  become  purified. 

One  commonly  used  method  for  obtaining  high  initial 
dilutions  is  to  have  a  series  of  ports  in  a  diffuser  pipe 
which  will  attach  to  the  end  of  the  outfall.  These  ports 
ensure  that  the  effluent  is  discharged  over  a  wide  area  thus 
ensuring  that  as  much  of  the  receiving  water  as  is 
practically  possible  will  be  used  to  dilute  and  render 
harmless  the  waste.  A  typical  diffuser  arrangement  is  shown 
in  Fig.  1-02.  In  certain  cases  diffusers  can  provide  initial 
dilutions  ten  or  more  times  greater  than  with  merely  an  open 
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TYPICAL  DIFFUSER  SECTION. 


FIG.  1-02.  COMMON 
DIFFUSER  ARRANGEMENT. 
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ended  outfall  ( Rawn  et  al,  I960).  Thus  nowadays  for  almost 
all  wastewater  disposal  schemes)  some  sort  of  diffuser  is 
recommended. 

Another  common  type  is  shown  in  Fig.  1—03.  It  consists 
of  a  buried  outfall  pipe  with  a  series  of  vertical  risers  on 
its  top  surface.  This  may  be  used  when  clearance  is  a 
problem,  for  instance  in  a  shallow  river,  as  the  main  body 
of  the  diffuser  is  buried  below  bed  level.  Any  excess 
velocity  that  an  effluent  Jet  may  initially  possess  will 
eventually  be  reduced  to  zero  and  if  it  has  buoyancy,  it 
will  then  tend  to  rise  to  the  surface.  In  the  case  of  a 
co— flowing  channel  where  the  effluent  and  the  receiving 
water  are  moving  in  the  same  direction,  a  plume  will  not 
usually  manifest  itself.  However,  the  effluent  may  still 
rise  to  form  a  surface  layer  for  some  distance  downstream  of 
the  point  of  injection. 

If  diffuser  ports  are  spaced  closely  together, 
interference  will  result  thus  diminishing  the  dilution 
capability  of  the  surrounding  water.  Almost  all  diffusers 
are  designed  with  the  aid  of  hydraulic  models  which  would 
indicate  whether  this  was  going  to  be  a  severe  problem.  It 
is  advantageous  to  have  equal  discharges  flowing  from  each 
nozzle  so  that  the  total  wastewater  load  may  be  spread  along 
the  whole  length  of  the  diffuser  pipe.  This  entails  careful 
design  of  both  the  exit  nozzles  and  the  pipe  itself.  One 
method  is  to  have  a  tapering  pipe  getting  smaller  towards 
its  end.  This  ensures  a  constant  velocity  inside  the  pipe 


, 


■ 

, 


7 


xz 


FIG.  1-03  DIFFUSER  WITH  RISER  PIPES 


FIG  1-04  BRIDGE  PIER  DISPOSAL 
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and  means  that  all  port  diameters  may  be  kept  the  same. 
Conversely  port  diameters  may  be  increased  towards  the  end 
of  the  pipe  for  the  same  effect.  Grace  (  1978  )  gives  a  method 
for  calculating  the  discharge  from  each  port  which  is  of  a 
step— by— step  nature  and  ideally  suited  to  computer 
appl ica t i on • 

Another  disposal  technique  Is  to  use  an  outfall  similar 
to  a  bridge— pier,  shown  in  Fig.  1-04,  with  outlets  on  either 
side  which  point  the  effluent  flow  across  the  path  of  the 
river  flow.  This  produces  two  distinct  trails  flowing  with 
the  river  which  will  eventually  meet  at  some  point 
downstre  am . 

A  modification  of  one  of  the  previous  systems  gives  an 
underground  pipe  with  risers  which  curve  through  forty— five 
degrees  to  lie  on  the  bed  of  the  channel  and  point 
downstream  thus  issuing  their  flow  in  the  same  direction  as 
the  receiving  water  flow.  Clay  and  Gordon  ( 1965 )  describe 
just  such  an  outfall  which  was  designed  for  a  pulp  and  paper 
mill  in  the  city  of  Kamloops,  B.C.  and  is  shown  in  Fig. 

1—05.  It  was  to  handle  a  discharge  of  0.65  cumecs  of 
effluent  to  the  Thomson  River,  a  salmon-producing  tributary 
of  the  Fraser  River.  Thus  it  was  extremely  important  that 
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DIRECTION  OF  RIVER  FLOW. 
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FIG.  1-05  CROSS-SECTION  THROUGH  THE  KAMLOOPS  DIFFUSER, 
SHOWING  A  TYPICAL  NOZZLE  ARRANGEMENT. 
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efficient  diffuser  through  which  it  was  passed.  This 
diffuser  provided  initial  dilutions  of  1:100  from  wh 
pollution  hazzard  would  arise. 

The  diffuser  itself  consisted  of  a  183  ra.  lengt 
cm.  diameter  pipe  lying  perpendicular  to  the  directi 
the  river  flow  buried  to  a  depth  of  30  cm.  To  this  p 
connected  110»  6.3  cm.  nozzles  each  joined  to  a  7.6 
diameter  rubber  hose,  2.45  m.  long.  Each  of  these  ho 
emits  a  3—D  jet  of  effluent  into  the  river  which  wil 
in  the  lateral  direction  towards  the  adjacent  jets  a 
in  the  vertical  direction  towards  the  water  surface. 

The  practical  work  for  this  thesis  consists  of 
into  the  velocity  distributions  at  different  vertica 
sections  downstream  of  one  nozzle  of  similar  design. 
Vertical  velocity  profiles  along  the  centre  of  the  j 
measured  as  were  lateral  profiles  which  indicate  eff 
spread  towards  adjacent  diffuser  ports.  In  a  practic 
at  some  distance  downstream  adjacent  flows  will  inte 
with  each  other  and  the  situation  will  be  reduced  to 
a  two  dimensional  slot  jet. 

Dilution  of  wastewater  that  occurs  downstream  o 
Injection  nozzle  is  a  direct  function  of  velocity 
distribution  and  analysis  of  jet  scales  shows  growth 
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II*  Theoretical  Considerations* 


Analysis  of  a  three-dimensional  jet  in  a  co— flowing 
channel  may  be  approached  by  first  considering  the  case  of  a 
circular  jet  with  an  efflux  velocity  of  Uo  in  an  infinite 
fluid  moving  with  a  velocity  of  Uj*  The  definition  sketch 
for  this  situation  is  shown  in  Fig*  2—01* 

This  diagram  shows  two  distinct  regions  of  flow*  Close 
to  the  nozzle  there  is  a  potential  core  or  flow  development 
region  followed  by  a  region  of  fully  developed  flow*  In  the 
fully  developed  flow  region,  the  velocity  profiles  are  shown 
and  these  are  found  to  be  of  a  similar  nature  at  every 
section  downstream  of  the  nozzle*  Scales  for  the  analysis  of 
the  profiles  are  now  introduced*  Um  is  the  excess  velocity 
of  the  jet  over  the  surrounding  flow  and  is  defined  by  Eq* 
2—01,  um  is  the  absolute  maximum  at  any  section,  Ui  is  the 
average  surrounding  flow  velocity,  the  x— axis  is  in  the 
downstream  direction  of  flow  and  the  r— axis  is  perpendicular 
to  this* 

Um  =  um— U i  2-01 

A  length  scale  b  is  now  introduced  where  b  is  the  distance 
from  the  point  of  maximum  velocity  to  a  point  where  the 
relative  velocity  is  half  the  maximum*  This  is  described  in 
Eq.  2-02. 

b  —  r(  Um)-r(  Um/2)  2-02 

The  equations  of  motion  and  the  continuity  equation  for  this 
jet  can  be  written  as  (Rajaratnam) 
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FID  201  DEFINITION  SKETCH  FOR  CIRCULAR  COMPOUND  JETS. 
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u9u  +  v9u_  _  1_  .  1_  9rx 
3x  3r  p  r  3r 


2-03 


3ru  3rv 

9x  3r 


2-04 


where  a  and  v  are  the  time  averaged  velocities  in  the  x  (or 
axial)  and  in  the  r  (or  radial)  directions  respectively,  p 
is  the  fluid  density  and  ^  is  the  turbulent  shear  stress* 

By  integrating  Eq*  2-03  with  respect  to  r  from  r= 0  to 
r=  oo  we  could  arrive  at  (  Ra  jaratnam,  1976  ) 

f°  purSu.dr  +  f°°  pvr9u.dr  _  1_  d_  /”  2irrdrpu2  1_  d_  f°°  2TTrdrpuUi 

0  3x  0  3r  2tt  dx  0  2tt  dx  0 

1_  d_  /°°  2iTrdrpu(u  -  Ui) 

2tt  dx  0 

r°°  . 

8rr.dr  r  -,«> 
and  J  o  =  {  rx  }  =0 

dr  0 

So  Eq*  2—03  becomes 


ld/°°  2xrdrpu(u  -  Ui)  _  n  2-05 

27d7  0  "  U 


which  says  that  the  rate  of  change  of  excess  momentum  flux 
in  the  x  or  downstream  direction  is  zero;  ie*  the  excess 
momentum  flux  of  the  jet  in  the  axial  direction  is 
preserved. 

The  similarity  between  velocity  profiles  at  different 
sections  downstream  of  the  nozzle  can  be  expressed  as 


Using  this  relation  along  with  Eq*  2—05  we  get 
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dx  0  ^U7 
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)dn 
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2-07 


In  the  case  of  a  jet  Issuing  into  a  fluid  that  is  travelling 
very  much  slower,  then  Um/Ui  close  to  the  nozzle  in  the 
fully  developed  flow  region  will  be  large.  This  situation 
where  the  jet  is  dominant  Is  known  as  the  strong  jet  case. 
Further  downstream  the  jet  effects  will  diminish  until  Urn  is 
of  the  same  order  as  (Jj  and  still  further  Urn  <<  Uj.  The 
latter  case  is  known  as  the  weak  jet  case. 

The  two  extreme  cases  could  now  be  considered 
separately  along  with  dimensional  analysis. 

For  the  strong  jet  where  Urn  >>  Eq.  2—07  becomes 

d  /°°b  2ri  (UiTh  2f2dn  =  0 

d7  0 

or 

d  (b2Um2)  /°°  T)2 f 2 dn  =  0  2-08 

dx  0 

The  integral  in  Eq •  2-08  will  be  a  constant,  therefore 

_ d(  b2Um2  )=0  2-09 

dx 

We  may  now  a  ssurae  that  both  "the  terms  within  th  e  brackets  in 
Eq •  2—09  are  some  function  of  x,  namely 

D  q 

Um  a  x  and  b  a  x 

Therefore  2p+2q=0  or  p+q=0  2—10 


For  the  weak  jet  case  Um  <<  U|  and  Eq.  2—07  becomes? 


- 
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d  /°°  b2r)Umfdri  =  0 

dx  0  Ux 


or  d  P  b2r|Umfdri  =  0 

d7  0 


Simplifying,  this  gives 


d  (b2Um)  f°°  riUmfdn  =  0 

dx  0 


and  _dX  b 2  U m  )  =  0  2-12 

dx 

and  therefore  2q+p=0  2-13 

When  Um  is  of  the  same  order  as  this  analysis  i  t 

applicable  and  other  methods  must  be  found. 

Now  using  the  similarity  of  the  equati  ons  o  f  Dio  t  i  on ,  we 
have  firstly; 

u  =  Ui  +  Um(  f  ) 


and 


u3u  =  Um'f  +  UmUm'f2  -  UiUmb 1  n f 1  -  Um2b_' nff' 
3x  b  b 


where  primes  denote  differentiation  with  respect  to  x. 
Integrating  the  continuity  equation  2—04  we  get 


and 

If 

then 


v  =  Umb ' 1  fnn2f 'dn 

o  0 


bUm' 1  /nofdn 
n  0 


9u 

3r 


Umf' 

b 


t/ (pUm2 )  =  q(n) 


3rx  _  5  (pUm2rg) 
3r  dr 


=  pUm2  (g  +  rg'/b) 


1_  1_  3rx 
p’r  3r 


2-14 


2-15 
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Using  these  expressions  along  with  Eq.  2—03  and  simplifying 


bUiUm'nf  .  bUm1  (nf2  -  f '/n  qfdn)  -  Ui_  b'n2f' 
Urn2  Um  0  Um 


-  b'(n2ff'  -  f 1  /n  n2f'dri)  =  g  +  ng '  2-16 

0 

If  we  now  assume  that  in  Ea.  2-16  all  terms  in  r|  ,  cr 
functions  thereof  are  the  same  order  of  magnitude,  then  the 
first  term  reduces  to, 


p 

q  xK 


1 


and  the  third  term  to 


UL  yP-1 
Um  x 


ie.  both  are  proportional  to  Utx 

Um 


q- 


For  the  strong  Jet  case  when  Um  >>  U*  the  second  and 
fourth  terms  in  the  equation  will  be  dominant, 


bUm'(nf2  -  f'/n  ufdn)  -  b'(n2ff'  -  f'/n  qf'dn)  =  a  +  ng ' 
Dm  0  0 


from  which  we  get 

bUm 1  a  x°  and  b*  a  x° • 

Dm 

This  gives  the  result  that  q-l=0  or  q=l.  Using  this  along 
with  Eq.  2—10  we  get 

q=l  and  p=— 1 


, 


- 


17 


This  means  that  for  a  strong  jet  Um  a  1/x 
and  b  at  x. 

For  the  weak  jet  the  second  and  fourth  terms  drop  out 
and  we  are  left  with 

bUiUm'nf  -  Uib'n2f'  =  g  +  nq1 
Urn2-  Um 


giving , 


bUm*  a  x°  and  b1  a  x° 

Um2  Um 

or  q— p- 1=0 

Using  this  result  along  with  Eq.  2—13  we  get 

q=l/3  and  p=— 2/3 

or  Um  at  1/x  2/3  an d  b  a  X1  / 3. 

We  have  already  demonstrated  that  excess  momentum  M0  is 
conserved  in  the  axial  direction.  Now  using  dimensional 
analysis  we  can  first  say  that 

Um  =  f  i  (  Mo  »p  »  x  ) 

wherein  the  viscous  effects  are  neglected,  which  is  valid  if 
the  jet  Reynolds  number  is  greater  than  a  few  thousand  and 
it  has  been  further  assumed  that  the  channel  is  wide  enough 
for  the  effects  of  the  side  walls  to  be  negligible  and  that 
the  flow  is  so  deep  that  the  presence  of  the  free  surface 
does  not  affect  t he  process  of  jet  diffusion.  Using  the 
Pi— theorem  we  get 


Um  Ci 

/(U0  ( U  o  -  U ! ) )  '  Wul 


Um 

/(M0/pxz ) 


Const. 


or 


- 
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or 

Um  Ci 

U,  '  iw 

is  the  so-called  momentum  "thickness  term  used  in  the 
of  scales  in  the  axial  direction  and  is  defined  by 
t  i  on  ; 

M0=  TT^O2.  pU0(  U0-U1  )=  T\e2.p\Jl2 

0  =  =  r0/(a(a-l ) ) 

is  the  initial  ratio  of  velocities  between  the  jet 
surrounding  fluid  (=Uq/Ui).  If  Uj  is  also  included 
nitial  list  of  independent  parameters  it  can  be 
at 

Um/U,  =  f  (x/6). 

Experimental  research  in  this  area  has  been  conducted 
by  Forstall  and  Shapiro  (  1950  ),  Tani  and  Kobashi  (1951)  and 
Antonia  and  Bi Iger  (  1973  )  among  others.  These  observers,  in 
trying  to  supply  a  formula  that  will  fit  their  data,  have 
chosen  a  variety  of  solutions.  Squire  and  Trouncer  ( 1944) 
proposed  that  a  relation  involving  the  cosine  of  (r/b)  be 
used  to  predict  velocity  in  the  fully  developed  region, 
while  Bradbury  and  others  indicate  a  power  relation  between 
all  the  relevant  terms.  Rajaratnam  (1976)  and  Pande  and 
Rajaratnam  (1976)  present  e  study  of  the  data  gathered  by 
some  of  the  observers  along  with  various  empirical  relations 


where  © 
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for  their  prediction#  While  each  observer  has  usually 
derived  a  relation  that  is  applicable  to  his  data, 
discrepancies  appear  when  all  the  data  is  correlated 
together  and  plotted.  Pande  and  Rajaratnam  have  developed  a 
method  of  analysis  using  both  the  integral  momentum  equation 
and  the  integral  energy  equation  that  predicts  U i /Um  for 
values  of  x/©  for  both  strong  and  weak  jets.  When  used  with 
the  results  by  Antonia  and  Bilger  this  equation  provides 
good  agreement.  It  is  also  used  in  the  analysis  in  Chapter  3 
of  this  thesis. 
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WALL  JETS. 
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the  wall  and  formerly  the  same  as  the  vertical 
increases  at  a  rate  approximately  four  to  five 
f  ast  • 

The  velocity  variation  in  the  inner  region 
approximated  by  a  power  law  distribution  of  the 


These  two  power  values  cover  the  extremes  of  an 
within  which  all  boundary  layer  results  would  g 
Fig*  3—12  shows  these  two  curves  drawn  on 
inner  region  velocity  variation  and  it  can  be  s 
almost  all  results  lie  inside  these  two* 

Ef  feet  Qf.  boun  da  rv  roughness 

During  all  the  data  collection*  the  flume 
artificial  roughness  on  the  bed  caused  by  a  str 
sandpaper,  which  is  described  in  more  detail  in 
Pani  (1972)  studied  the  effect  of  boundary  roug 
wall  Jet  in  a  stagnant  receiving  environment*  K 
measurements  obtained  for  jet  velocities  over  a 
aluminum  bed  with  those  for  the  same  nozzle  ove 
identical  sandpaper  to  that  used  in  this  resear 
He  concluded  that  the  presence  of  a  rough 
not  affect  the  maximum  velocity  decay  along  the 
of  the  Jet  to  any  appreciable  extent*  Plotting 
smooth  wall  data  together  in  a  Umo/U0  vs*  x//  A 
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d/ks=  1,  where  ks  is  the  equivalent  sand  roughne 
sandpaper*  He  considered  two  distinct  regions  of 
first  being  the  inner  or  boundary  layer  where  th 
increases  to  its  maximum  value  and  the  second  be 
outer  or  free  mixing  region  where  the  velocity  d 
maximum  to  zero* 

Analysis  showed  that  the  presence  of  roughn 
the  boundary  layer  to  increase  In  size  and  thus 
velocity  to  occur  further  from  the  wall*  However 
mixing  region  remains  unaffected  by  the  presence 
roughness*  Pani  also  found  that  the  roughness  ha 
effect  on  the  inner  region  close  to  the  nozzle* 
be  expected  as  here  the  boundary  layer  is  thin  a 
of  roughness  height  to  boundary  layer  size  is  mu 
than  for  sections  further  downstream*  He  conclud 
presence  of  boundary  roughness  has  no  effect  on 
spread  of  the  jet* 

As  most  real  situations  such  as  dicharge  th 
diffuser  into  a  river  include  the  presence  of  bo 
roughness  and  as,  according  to  the  work  of  Pani, 
roughness  did  not  significantly  affect  the  analy 
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Ca  Icula  t  ion  of  dilution  from  vel  oc  i  tv  dlstri  but  j.  on 

If  we  consider  the  case  of  a  pollutant  added  to  a 
moving  flow  then  its  concentration  will  diminish  as  it 
spreads  out  in  all  directions  downstream  due  to  turbulent 
mixing,  entrainment  and  ambient  turbulence# 

An  analysis  similar  to  that  performed  on  the  jet  in  the 
first  section  yields  the  result  that  for  a  3— D  pollutant 
source,  pollutant  concentration  at  any  point  can  be 
expressed  in  terms  of  the  maximum  along  the  diffuser 
centreline  and  its  distance  from  the  centreline*  This 
situation  is  entirely  analagous  to  the  spreading  of  the  jet 
and  we  can  say  that 


c_ 

Cm 

=  f(n) 

and 

Cm 

s 

a  x  , 

b  a  x^  , 

where  Cm  is 

the 

maximum 

concentration  at  any  section 

downstream  of  the  source,  C  is  the  concentration  at  any 
point,  z  is  the  distance  from  the  centreline  and  b  is  the 
distance  to  half  the  maximum  concentration* 

We  may  again  say  that  Cm  and  b  are  both  functions  of  x 
and  using  the  pollutant  conservation  equation  we  can  show 
that  the  variation  of  Cm  in  the  downstream  direction  is  the 
same  function  of  x  as  the  velocity  scale  Urn.  From  these 
considerations,  we  could  Indicate  that  for  the  strong  as 
well  as  weak  jet  cases  the  concentration  will  vary  in  the 
same  manner  as  the  velocity  scale*  Of  course  these 


proportionalities  do  not  contain  any  constants  and  this  is 
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where  the  relations  will  differ#  This  is  comparable  to  the 
case  of  a  passive  plume  where  a  Jet  is  injected  into  a 
channel  with  zero  excess  momentum#  The  equation 

f  (n)  =  e-°'693rl2 


can  be  used  to  describe  the  profiles  of  a  pollutant  in  both 
the  outer  region  and  the  transverse  direction,  spreading 
downstream  of  a  source# 

It  has  been  found  by  experimentation  that  the 
concentration  spreading  from  a  circular  turbulent  jet  is  15 
to  20  percent  greater  than  the  velocity  spread.  That  is  the 
length  scale  b  is  greater  by  that  amount  at  every  downstream 
section  for  the  pollutant#  Various  researchers  have 
conducted  experiments  into  the  exact  value  of  this 
difference  and  among  them  the  value  of  17  percent  by  Hinze  & 
Zljnen(1949)  appears  to  be  the  best. 

If  we  include  this  factor  in  the  equation  it  becomes 
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' 


25 


So  the  rate  at  which  a  pollutant  is  diluted  can  be 
determined  for  a  particular  situation  and  from  the 
guidelines  that  have  been  proposed  for  effluent  disposal  it 
may  also  be  determined  whether  the  propsed  installation  will 
prevent  pollution  of  the  receiving  water. 


oJL  idom£JO.tilQi  XliiX  J.H  the  downs  tream  dj  reel  ion 
In  the  original  analysis  it  was  shown  that  the  momentum 
efflux  in  the  downstream  direction  was  preserved.  This  is 
not,  in  fact,  strictly  true  as  the  shear  stress  on  the  bed 
will  cause  the  efflux  to  diminish. 

If  the  momentum  efflux  is  decreased  by  an  appreciable 
amount,  say  more  than  25  oercent,  then  the  analysis  is  no 
longer  applicable  to  sections  downstream  of  this  area.  The 
actual  decrease  was  therefore  computed  for  runs  5  and  7 
which  have  two  different  nozzle  diameters. 

The  procedure  used  was  an  approximate  one  but  was 
assumed  to  be  reasonably  accurate.  Instead  of  an  exponential 
variation  of  velocity  in  the  lateral  direction,  a  straight 
line  was  assumed  and  although  it  has  been  shown  that  the 
actual  spread  of  the  jet  is  about  2.25bz  a  value  of  2bz  was 
used  as  the  jet  is  an  elliptical  shape  and  its  width  would 
tend  to  diminish  as  it  approached  the  bed. 

Drag  force  was  computed  by  splitting  the  downstream 
area  into  twelve  sections.  For  each  section  an  average  shear 
stress  was  calculated  along  the  centreline  and  this  was 


again  averaged  along  with  the  free  stream  shear  stress.  The 
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area  of  the  trapezoidal  section  on  one  side  of  t 
determined  and  the  two  values  multiplied  togethe 
the  average  shear  force  for  one  side  of  a  sectio 

The  drag  force  at  each  other  section  was  fo 
these  were  summed  together  and  doubled  to  give  t 
drag  force  on  both  sides  of  the  jet.  The  percent 
total  momentum  excess  was  then  determined  at  the 
sect  ion. 

It  was  found  that  the  figure  of  25  percent 
in  both  cases  at  the  last  downstream  section.  Wi 
nozzle  and  a  =  19.89f  the  momentum  reduction  was 
at  the  last  section  and  for  run  7  with  the  large 
a  =  3.54  it  was  40  percent. 

However  for  the  first  case  the  reduction  wa 
limits  up  to  the  eleventh  section  and  for  the  se 
was  within  limits  up  to  the  tenth  section. 

It  was  therefore  felt  that  although  a  reduc 
it  was  not  significantly  large  in  the  main  resea 
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Ill#  Experimental  Arrangement  and  Experimental  Results# 


The  apparatus  used  in  this  study  consisted  of  a  flume 
approximately  50  metres  in  length  and  just  under  one  metre 
in  width#  The  slope  of  the  flume  can  be  varied  by  use  of  a 
motor  and  screw  jacks  at  10  metre  intervals  but  for  all 
experiments  the  slope  was  kept  at  a  constant  value  of  0#03%# 
A  tail  gate  was  used  to  control  water  level  in  the  flume  and 
a  valve  on  the  water  supply  line  altered  the  discharge  from 
0  to  2«5  c  «  f • s • 

Along  the  bed  of  the  flume  for  a  distance  of  15  metres 
cloth  sandpaper  was  glued  to  give  an  artificial  roughness# 
This  sandpaper  is  described  in  greater  detail  in  the  section 
on  shear  stress# 

Inject  i.on  JQZ^JLe 

The  method  chosen  for  water  Injection  was  by  two 
different  diameters  of  nozzles#  The  first  had  an  exit 
diameter  of  1/4  inch  and  was  used  for  Runs  1  through  6  and  a 
large  nozzle  of  3/4  inch  diameter  was  used  for  Runs  7  and  8# 
The  alternatives  in  the  placement  of  the  nozzle  consisted  of 
either  cutting  a  hole  in  the  bottom  of  the  flume  and  have 
the  water  supply  come  in  through  the  bed  or  else  suspend  the 
nozzle  from  some  point  above  the  bed  in  order  to  avoid 
drilling  holes  in  the  flume#  The  latter  alernative  was 
rejected  as  it  was  felt  that  a  supply  pipe  to  the  nozzle 
coming  down  through  the  flow  would  seriously  upset 
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downstream  velocity  distributions  by  causing  turbulent 
eddies  behind  the  pipe*  It  was  therefore  necessary  to  cut  a 
hole  2*5  inches  in  diameter  in  the  bottom  of  the  flume  to 
accomodate  the  nozzle  as  shown  in  Fig*  3—01* 

In  order  to  minimize  flow  disturbance  as  the  water  is 
pointed  downstream)  the  angle  of  the  elbow  upstream  of  the 
nozzle  was  made  45  degrees*  This  entailed  a  larger  cut  out 
in  the  flume  bed  but  was  considered  to  be  worth  the  extra 
trouble  to  avoid  irregular  velocity  patterns  in  the  flow 
created  by  a  right  angled  bend* 

£oflstaat  bead  arrao^ejent 

A  simple  constant  head  tank  was  constructed  to  maintain 
a  steady  flow  of  water  into  the  nozzle  at  a  constant 
velocity*  This  tank  was  mounted  on  the  side  of  the  flume  as 
shown  in  Fig*  3—02  and  had  a  fine  adjustment  that  enabled 
accurate  heads  to  be  maintained*  The  head  difference  between 
the  water  level  in  the  flume  and  that  in  the  constant  head 
tank  was  measured  by  a  point  guage  with  a  vernier  scale* 

Originally  the  water  supply  to  the  nozzle  came  direct 
from  the  city  fresh  water  supply*  This  was  found  to  be 
satisfactory  for  the  first  and  third  runs  but 
inconsistencies  were  observed  in  the  second  run*  It  was 
found,  while  measuring  velocities  at  various  sections 
downstream  of  the  nozzle,  that  the  jet  would  apparently 
disappear  and  just  as  suddenly  then  reappear* 


The  cause  of  this  problem  remained  a  mystery  for  a  long 
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Direction  of  Flow.  - ► 


INJECTION  NOZZLE. 


Jet  water 
from  sump. 


Overflow. 


FIG.  3-02.  CONSTANT  HEAD  ARRANGEMENT. 
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time  until  it  was  noticed  that  an  air  lock  would  appear  in 
the  nozzle  supply  line  on  occasions.  Further  investigation 
revealed  that  as  the  city  water  came  into  the  the  warm 
laboratory,  air  bubbles  would  form  and  these  would 
eventually  Join  together  to  form  one  large  air  lock  at  the 
beginning  of  the  nozzle  suoply  line  thus  interrupting  flow 
to  the  jet.  The  reason  for  this  problem  not  arising  in  the 
other  runs  was  that,  in  the  first  case,  the  jet  disc  barge 
was  sufficiently  high  to  carry  air  bubbles  into  the  flume 
through  the  nozzle  while  in  the  third  run  the  discharge  was 
sufficiently  low  to  allow  the  bubbles  to  rise  to  the  surface 
under  their  own  buoyancy. 

After  the  cause  was  discovered,  the  city  supply  was 
discontinued  and  instead  a  submersible  pump  was  placed  in 
the  water  sump  to  supply  the  constant  head  tank.  This 
appears  to  be  a  better  solution  as  temperature  and  density 
differences  between  the  flume  water  and  jet  water  would 
cause  only  negligible  buoyancy  effects.  No  further  problems 
with  the  appearance  of  air  bubbles  were  experienced  as  the 
sump  water  had  become  de— aired  over  a  long  period  of  time. 


, 
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Velocity  Mep. guremgn 

All  velocity  measurements  were  made  with  a  pitot— static 
tube  1/16  inch  in  diameter.  It  was  considered  that  this  tube 
was  of  small  enough  size  so  as  not  to  cause  severe 
disturbance  to  the  flow.  The  two  outlets  from  the  pitot  tube 
were  connected  to  either  side  of  a  pressure  transducer  which 
sensed  the  magnitude  of  the  differential  pre  ssu  re  at  any 
point  in  the  flow.  This  pressure  is  entirely  due  to  the 
velocity  at  that  point. 

The  transducer  consists  of  a  magnetic  pressure  plate 
which  is  moved  towards  one  side  of  the  transducer  body  by 
the  exc  ess  head  co  rrespondi ng  to  the  velocity  at  a  certain 
point.  This  magnetic  signal  Is  interpreted  by  a  transducer 
meter  and  output  as  a  voltage. 

The  calibration  of  the  transducer  was  accomplished  by 
two  water  supply  pots.  By  the  use  of  a  valve  board,  the 
pressure  supply  can  be  switched  between  the  calibration  pots 
and  the  pitot  tube  itself.  The  transducer  is  first  zeroes  oy 
having  identical  heads  in  both  calibration  pots  and  then  a 
differential  head  of  1  inch  is  aoplied  to  the  transducer  and 
the  meter  adjusted  till  this  corresponds  to  an  output  of  one 
volt.  Thus  a  voltmeter  would  indicate  velocity  head  directly 
in  inches  of  water. 

It  is  extremely  important  that  no  traces  of  air  remain 
in  the  lines.  For  this  purpose  de— aired  water  was  used.  Any 
air  bubbles  that  exist  in  the  lines  cause  n on-l i nea r 1 ty  of 
head  measurements  as  the  air  bubbles  become  compressed. 
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After  calibration*  the  supply  pots  are  switched  out  and  the 


transducer  is  connected  to  the  pitot  tube  which  has  also 
been  carefully  de— aired# 

sh e^r  Strega  avajuaijon 

The  measurement  of  dynamic  head  on  the  bed  of  the  flume 
at  various  distances  downstream  of  the  injection  nozzle  was 
made  at  the  time  of  the  general  velocity  distribution 
measurement  and  from  this  head  the  shear  stress  was  later 
calculat  ed • 

This  method  was  suggested  by  Preston  in  1954#  Preston 
assumed  that  the  total  pressure  P  indicated  by  a  total  head 
tube  resting  on  the  wall  and  facing  the  flow  minus  the 
static  pressure  at  the  wall  Po  is  a  function  of  only  the 
shear  stress  at  the  wall  Tot  diameter  of  the  tube  d,  the 
mass  density  p  and  the  kinematic  viscosity  of  the  fluid  \) , 
Using  the  Pi— theorem,  the  above  statement  could  be  reduced 
to  , 


Patel  ( 1965)  found  Preston's  original  equations  to  be 


somewhat  in  error  and  suggested  three  new  relations 


connect i ng 


Ap*  =  Apd2 
4pv^ 


3-02 


and 


, 


- 


TO*  -  Tpd2 

4pv2 

for  smooth  boundaries*  Hwang  and  Laursen  (1963)  appl 
Preston’s  work  to  rough  boundaries  and  found  that,  w 
modification,  the  technique  was  applicable  there  als 
In  the  case  of  rough  boundary  flow,  further  par 
must  be  introduced,  these  being  u* ,  the  friction  vel 
and  k  the  roughness  height  of  the  bed.  For  completen 
internal  diameter  of  the  tube  opening,  a,  could  be 
introduced  but  it  was  demonstrated  by  Hsu  (  1955  )  tha 
effect  of  including  this  parameter  was  small* 

Again  using  dimensional  analysis  these  groups  c 
obta ined , 

Ap  =  f  ,u*k  ,  a) 

^  U  V 

o  r  Ad  =  f  |U*k  ,  d 

to  1  v  kJ 

In  the  fully  rough  flow  region,  viscosity  is  no  long 
Important  and  the  ratio  of  dynamic  pressure  to  shear 
becomes  a  function  of  tube  diameter  to  roughness  siz 
Furthermore,  if  one  tube  is  used  along  with  constant 
roughness,  d/k  itself  becomes  a  constant* 

Hollingshead  (1972)  produced  a  graph  of 


where  ks  is  Nikuradse*  s  equivalent  sand  grain  roughn 
using  the  equations  of  Patel  and  Nikuradse  and  his 
experimental  measurements*  This  curve  is  shown  in  Fi 
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PRESTON  TUBE  CALIBRATION  ANALYTICAL  SOLUTION 
WITH  OTHER  EXPERIMENTAL  DATA 


FIG.  3-03 


GRAPH  FOR  SHEAR  STRESS 


PREDICTION. 
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The  roughness  used  throughout  the  testing  was  a  length 
of  No*  36  aluminum  oxide  wet  or  dry  cloth  shown  in  Plate 
3—01*  It  was  placed  on  the  bed  of  the  flume  at  a  sufficient 
distance  upstream  of  the  measurement  area  to  insure  a  fully 
developed  rough  wall  boundary  layer.  Hollingshead  (1972) 
found  that  the  grain  size  was  0.00183  feet  which  gave  an 
equivalent  roughness  ks  of  0.00607  feet.  The  curve  used  for 
shear  stress  prediction  was  the  one  for  a  d/ks  value  of 
0.0067/0.00607  or  approximately  1. 

Hollingshead's  graph  gives  a  very  complete  treatment  of 
shear  stress  determination  in  the  smooth  turbulent  flow, 
transition  and  rough  turbulent  flow  regions  and  provides  a 
relatively  simple  method  of  finding  the  wall  shear  stress 
from  easily  raeasureable  parameters.  All  shear  stress  points 
for  this  research  lay  on  the  constant  part  of  the  curve  for 
fully  developed  rough  turbulent  flow. 

in^tr.^aeniaAlQD 

Data  collection  was  by  means  of  a  system  originally 
developed  to  study  3— d i mens  1 ona 1  flow  around  a  cantilevered 
obstruction.  It  consisted  of  a  control  unit,  traverse  unit 
and  logging  and  recording  unit.  The  control  unit  was  a 
microprocessor  activated  3— D  traverse  positioner  containing 
a  computer  program  written  in  basic  language.  This  program 
allowed  various  modes  of  operation  from  numerical  point 
entry  to  coordinate  editing  to  locating  the  bed  of  the  flume 


at  various  measurement  sections.  All  Information  entry  was 
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All  divisions  are  in  millimetres. 
Plate  3-01.  BED  ROUGHNESS. 
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via  a  teletype  and  usually  in  the  form  of  punched  pa 
for  hard  copy  storage. 

A  typical  run  necessitated  firstly  loading  the 
program  depending  on  whether  a  lateral  or  vertical  s 
to  be  performed,  then  entering  the  coordinates  that 
pitot  tube  was  to  visit.  After  entering  the  run  comm 
microprocessor  would  issue  signals  to  a  motorized  tr 
whose  exact  location  was  measured  by  potentiometers. 

The  traverse  unit  contains  a  carriage  for  the  p 
tube  mounted  on  stainless  steel  rods  and  using  linea 
bearings  for  motion  in  the  vertical  and  lateral  dire 
Potentio  me  ters  on  the  non— slip  driving  cables  e n  sure 
location  of  the  tube  in  all  directions  every  time  th 
measurements  are  repeated.  In  the  downstream  di recti 
traverse  unit  rides  on  three  aluminum  wheels  and  tra 
the  top  of  the  flume. 

Once  the  pitot  tube  and  traverse  unit  have  been 
positioned  at  a  particular  co-ordinate  a  delay  of 
approximately  forty— five  seconds  ensues  in  order  to 
conditions  in  and  around  the  pitot  tube  time  to  sett 

The  coordinates  and  transducer  values  are  scann 
Fluke  2240A  datalogger  with  60  channels  entered.  The 
four  channels  read  are  the  y,  x,  and  z  coordinates  a 
water  level  reading  and  the  remaining  56  channels  ar 
repeated  readings  of  the  pressure  transducer  for  the 
tube  with  periodic  readings  of  the  water  level  in  be 

The  whole  process  of  datalogging  takes  about  35 
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for  each  coordinate,  so  it  can  be  seen  that  an  avera 
value  of  the  transducer  readings  will  yield  a  very  a 


result.  All  60  channels  are  recorded  on  a  digital 


ca 


for  each  coordinate  and  once  the  run  is  started,  the 
traverse  unit  will  keep  moving  from  point  to  point, 
its  data  each  time,  till  either  the  cassette  is  full 
coordinates  have  been  visited. 

The  traverse  and  control  equipment  are  illustra 
Plate  3—02.  Once  the  cassette  has  been  filled,  it  is 
transferred  to  a  PDP/8  comouter  which  has  a  program, 
in  basic,  to  reduce  the  60  channels  to  one  line  of  6 
readings.  Those  6  readings  are  the  y,  x,  and  z  coord 
the  water  level,  the  transducer  pressure,  averaged  f 
readings,  and  the  corresponding  velocity  at  that  poi 

Once  again  a  typical  run  consists  of  12  section 
downstream  of  the  nozzle  with  usually  21  measured 
coordinates  in  either  the  vertical  or  horizontal  dir 
so  altogether  454  points  were  visited  for  each  run  a 
velocity  measured  at  all  of  them. 

All  the  data  from  the  PDP/8  was  then  transferre 


main  university  system  for  subsequent  plotting  and  f 
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Plate  3-02.  Traverse  and  data  aquisition  system. 
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Experiment jej:og.e^ar.e  results 

The  practical  work  for  this  thesis  was  mainly  composed 
of  eight  separate  series  of  velocity  distribution 
measurements*  The  first  six  experiments  were  performed  with 
Jet  injection  through  a  0*25  inch  nozzle  and  the  last  two 
through  a  nozzle  of  0*75  inch  diameter. 

Various  ratios  of  velocity  between  the  jet  and  the 
surrounding  channel  flow  were  chosen  to  cover  as  broad  a 
spectrum  as  possible  and  these  are  listed  in  Table  3—01*  A 
suitable  channel  discharge  was  selected  and  with  the  flume 
width  constant  at  2*97  feet,  a  convenient  depth  was  chosen* 
This  depth  was  to  be  kept  at  0*61  feet  throughout  the 
research  as  it  was  felt  that  this  gave  an  adequate  depth  to 
nozzle  size  ratio  so  as  not  to  cause  any  free  surface 
effects  on  the  jet* 

The  constant  head  tank  was  then  adjusted  to  provide  the 
chosen  velocity  ratio*  The  calculation  of  jet  velocity  was 
first  taken  as  /2gh  where  h  is  the  dlf  ference  i n  he i gh  t 
between  the  wa  te  r  surface  i n  the  constant  head  tank  and  the 
free  stream  surface  and  g  is  the  acceleration  due  to 
gravity*  On  measuring  the  actual  velocity  of  flow  from  the 
nozzle  with  the  pitot  tube  however,  it  was  found  to  be 
considerably  lower,  due  to  the  friction  and  fitting  losses 
in  the  nozzle  feed  pipe*  Therefore  the  actual  velocity  ratio 
was  calculated  using  the  measured  velocity  of  the  jet  rather 
than  the  calculated  value* 


Sections  at  suitable  distances  downstream  of  the  nozzle 
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were  chosen  for  velocity  distribution  measurement  and  these 
are  shown  in  the  graphs  in  Appendix  B  In  an  x/d  form#  The 
microprocessor  deals  not  in  absolute  distances  but  rather  in 
co-ordinates  which  translate  into  voltage  measurements  at 
the  potentiometers  located  on  the  traverse  lines*  Thus  the 
chosen  distances  were  converted  Into  their  respective 
co-ordinates  before  entry  into  the  computer.  After  selection 
of  the  points  of  velocity  measurement,  the  next  procedure 
was  to  locate  the  exact  bed  of  the  flume  at  these  points. 
This  is  described  in  the  section  on  Instrumentation. 

After  a  co-ordinate  corresponding  to  the  bed  position 
was  found,  a  choice  had  to  be  made  of  how  many  positions  in 
the  vertical  direction  were  to  be  measured.  For  the  first 
run  a  figure  of  38  was  chosen  and  the  values  of  the 
co-ordinates  that  would  be  added  to  the  bed  value  worked 
out.  Because  of  the  large  velocity  gradient  in  the  lower 
region  of  the  flow,  the  co-ordinates  were  naturally  spaced 
more  closely  together  near  the  bed. 

After  all  values  were  selected,  the  run  command  was 
entered  and  the  microprocessor  took  over.  It  first  issued 
commands  to  move  the  traverse  to  the  section  closest  to  the 
nozzle  and  along  its  centreline,  then  dropped  to  the  bed  and 
after  a  suitable  delay  began  scanning  the  data  and  recording 
it  on  cassette  tape.  On  completion  of  measurement  at  the 
bed,  it  moved  to  another  section  above  it  and  repeated  the 
process.  Because  some  of  the  vertical  points  of  measurement 
were  so  close  together,  the  movement  of  the  pitot  tube 


, 
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carriage  was  in  an  up  and  down  fashion  enabling  it  to  move 
and  come  to  rest  accurately  at  the  proper  location. 

It  was  found  that  38  points  in  the  vertical  direction 
were  more  than  enough  to  obtain  an  accurate  profile  of 
velocity  variation  and  this  was  subsequently  reduced  to  27 
points  and  then  again  to  21  points. 

Cnee  the  pitot  tube  had  visited  all  vertical  points  at 
one  section,  it  then  moved  downstream  to  the  next  one,  while 
still  remaining  on  the  nozzle  centreline  and  repeated  the 
whole  process.  After  6  sections  the  cassette  tape  was 
usually  full  and  had  to  be  changed.  The  data  on  the  tape  was 
then  analysed  and  reduced  as  described  previously. 

From  this,  the  point  of  maximum  velocity  occurence 
could  be  easily  found  and  this  was  the  vertical  position 
chosen  for  the  lateral  scan.  The  program  in  the 
microprocessor  was  modified  slightly  to  scan  in  this  new 
direction  and  suitable  incremental  co-ordinates  were  again 
chosen. 

The  lateral  scan  was  started  slightly  to  the  left  of 
centre,  to  check  on  symmetry  of  the  Jet  as  most  measurements 
were  confined  to  one  side  alone. 

The  procedure  was  repeated  for  the  same  12  sections 
downstream  of  the  nozzle  and  after  all  data  had  been 
collected  it  was  known  accurately  how  the  Jet  decayed  along 
the  centreline  of  the  nozzle  and  how  it  spread  in  the 


lateral  direction 
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Ef  feet  o  f  de p th  a  l tera  t i on 

Although  the  depth  of  flow  was  kept  constant  throughout 
the  experiments  two  different  depth  to  nozzle  size  ratios 
were  investigated#  It  was  decided  that  a  study  of  the  effect 
of  flow  depth  alteration  on  the  Injected  water  would 
complete  the  research#  Accordingly  a  standard  velocity  ratio 
of  10:1  was  chosen  for  all  tests  and  the  depth  was  reduced 
step  by  step  while  photographs  were  taken  for  later 
comparison#  When  these  photographs  were  studied  together,  it 
was  noticed  that  little  or  no  effect  was  exerted  on  the 
plume  by  the  depth  alteration#  The  actual  reduction  was  from 
the  usual  depth  of  0*61  ft#  down  to  0#22  ft#  which 
corresponds  to  a  depth  to  nozzle  ratio  of  3*25:1#  At  the 
lowest  depth  considerable  rippling  of  the  water  about  1  ft# 
downstream  of  the  nozzle  was  noticed  indicating  that  the  jet 
water  had  reached  the  surface#  During  the  experiments 
conducted  at  the  normal  water  depth  the  position  where  the 
jet  water  finally  surfaced  was  not  noticed  till  a  large 
distance  downstream#  The  scales  layed  on  the  bed  of  the 
flume  Indicated  the  lateral  spread  to  be  approximately  the 
same  for  all  depths.  The  depth  to  nozzle  size  ratio  varied 
from  about  10:1  down  to  3:1.  It  is  safe  to  say  that  there 
will  be  few  installations  where  the  above-mentioned  ratio 
will  approach  or  go  below  10:1# 

It  should  be  mentioned,  however,  that  this 
investigation  was  confined  to  a  short  distance  downstream  of 


the  nozzle  and  with  deoth  reduction  the  plume  will  naturally 
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reach  the  the  water  surface  much  sooner  at  some  point 
downs  tre  am • 

Aga.lysja  el  gyjsejrj. mental  r.s.su.ljt 

Most  analysis  and  reduction  of  the  data  was  done  on  the 
main  University  computer.  A  Listing  of  aLL  co-ordinates  and 
velocities  for  each  point  visited  by  the  pitot  tube  was 
obtained  and  from  this  the  maximum  velocity  ( um )  for  each 
section  was  read,  the  height  above  the  bed  at  which  this 
occured  was  found  and  the  surrounding  velocity  of  flow 
calculated  from  averaging  the  Last  6  or  7  readings  near  the 
surface.  This  average  velocity  was  used  in  calculations 
rather  than  the  one  read  from  the  orifice  meter  in  case  of 
any  slight  cal i bra  tion  errors  in  the  t  ransducer • 

The  co-ordinates  first  obtained  from  the  datalogger 
were  only  voltages  and  so  they  had  to  be  converted  to  a 
suitable  length  unit,  in  this  case  feet.  Next  the  value  of 
half  the  relative  velocity  was  calculated  and  its  point  of 
occurence  found  by  Interpolation  between  two  positions.  Once 
this  point  was  found  the  value  of  6  was  subtracted  from  it 
and  this  then  became  the  b  value  for  length  scale  analysis. 

This  procedure  was  repeated  for  every  downstream 
section  and  thus  the  four  parameters  used  in  the  analysis, 
the  maximum  velocity  um,  the  surrounding  flow  velocity  Ut, 
the  height  of  maximum  velocity  occurence  6  and  the  distance 
to  half  the  maximum  velocity  b,  were  obtained. 

For  the  lateral  spread  the  analysis  was  similar  except 
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illustration  here.  The  rest  of  these  plots  are 
Appendix  B.  Dimensional  plots  of  the  twelve  ver 
lateral  sections  are  shown  in  Fig's.  3  —  04  to  3—0 
The  variation  of  the  length  scales  in  the 
direction  is  now  illustrated.  Fig.  3—08  shows  t 
the  boundary  layer  thickness  and  it  can  be  seen 
through  6,  with  the  small  nozzle  all  lie  in  app 
the  same  area  while  the  values  for  the  larger  n 
to  decrease  in  the  downstream  direction.  This  i 
to  the  nozzle  being  placed  in  a  slightly  downwa 
affecting  the  normal  increase  in  boundary  layer 
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FIG.  3-10  Growth  of  transverse  length  scale  in  the  axial  direction 
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ALPHA  =  15:1  ALPHA  =  20:1  ALPHA  =  34:1 

Plate  3-03  Growth  of  jet  for  different  velocity  scales. 
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the  outer  or  free  mixing  region,  the  Inner  or  boundary  layer 
region  and  the  transverse  region#  For  the  outer  region  a 
non-dimensional  plot  of  (  y^-  §)/b  vs*  U/Um  was  produced,  where 
the  velocities  are  with  respect  to  the  free  stream*  Also 
plotted  was  the  exponential  curve  and  it  was  found  that  very 
close  agreement  exists  between  the  data  from  this  research 
and  this  curve*  All  data  in  the  boundary  layer  and  more  than 
a  ( y— 6)/b  value  of  about  3  was  ignored  in  this  plot  which  is 
shown  in  Fig*  3—11* 

In  the  boundary  layer  y/6  was  plotted  against  u/um 
where  in  this  case  the  velocities  are  absolute.  The  data  is 
found  to  lie  within  a  power  law  envelope  as  described 
earlier  and  shown  in  Fig*  3—12* 

In  the  case  of  Runs  1  through  6  the  first  two  sections 
downstream  of  the  nozzle  were  not  included  as  these  lay 
within  the  potential  core  region.  When  plotted  up  along  with 
the  other  sections  they  displayed  considerable  inconsistency 
as  the  Jet  was  not  developed.  In  the  runs  of  lower  velocity 
ratio  the  jet  effects  had  diminished  a  short  distance 
downstream  and  some  of  the  further  sections  were  omitted* 

The  velocity  variation  in  the  lateral  direction 
although  different  from  that  in  the  upper  region  can  again 
be  described  by  an  exponential  curve  at  each  section*  The 
lateral  spread  from  the  data  of  Run  S  is  shown  in  Fig.  3—13* 
The  remainder  of  these  plots  for  the  other  velocity  ratios 
appear  in  Appendix  B  also* 


The  scales  used  in  non— dl me ns i onal i si ng  the  previous 
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plots  were  then  analysed  along  with  other  terms  such  as  the 
momentum  thickness  0»  the  shear  velocity  U*  and  the  bed 
shear  stress  t  o« 

Shown  in  Fig.  3—14  is  the  non-dimensional  variation  of 
maximum  centreline  velocity  Urn  in  the  x— direction.  The 
momentum  thickness  term  0  is  used  as  the  axial  length  scale 
and  this  then  takes  into  account  any  large  momentum  effluxes 
so  that  data  from  large  or  small  velocity  ratios  will  plot 
along  the  same  line.  It  can  be  seen  that  agreement  is  good 
in  the  lower  areas  and  ignoring  the  data  of  Run  3,  a  mean 
line  can  be  drawn.  The  results  of  Run  3  have  caused  problems 
before  due  to  the  very  low  velocity  ratio  and  the  difficulty 
of  detecting  any  velocity  variation  far  from  the  nozzle.  It 
is  therefore  safe  to  say  that  the  jet  effects  in  this  run 
have  dissappeared  a  short  distance  downstream  of  the  nozzle 
and  it  is  valid  to  ignore  points  further  away.  The  other  two 
curves  shown  on  this  graph  are  the  results  of  Rajaratnam  for 
similar  analysis  applied  to  plane  wall  jets  and  circular 
jetsf  both  in  co— flowing  channels.  The  agreement  is  good  up 
to  an  x/0  value  of  about  seven  and  after  that  the  effect  of 
the  wall  becomes  significant  and  the  maximum  velocity  Urn 
decreases  much  more  quickly. 

Fig.  3-15  shows  the  variation  of  bed  shear  stress  t  o  -n 
the  axial  direction.  Some  of  the  measurements  were  made  very- 
close  to  the  nozzle  before  the  jet  had  time  to  properly 
establish  itself  on  the  bed  against  the  wake  effects  of  the 


nozzle  body  and  this  is  the  reason  for  the  low  shear  stress 


' 


' 


Expt.  1  12.62 

2  4.47 

3  2.21 

4  7.03 


60 


Prediction  of  velocity  scale  at  a  downstream  centreline  section 
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x(in) 

FIG.  3-15  Variation  of  bed  shear  stress  in  the  axial  direction 
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values  near  -the  nozzle# 

Figs#  3—16  and  3—17  show  the  non-dimensional  variation 
of  the  vertical  and  lateral  length  scales.  It  can  be  seen 
that  the  lateral  scale  grows  significantly  more  slowly  than 
its  counterpart  in  a  stagnant  receiving  environment#  This  is 
to  be  expected  as  the  flow  in  the  channel  past  the  nozzle 
restricts  its  growth# 

Fig#  3—18  describes  the  variation  of  shear  velocity  U# 
in  the  axial  direction#  Good  agreement  is  found  between  the 
data  for  all  runs  but  there  Is  a  significant  difference 
between  this  curve  and  that  for  plane  wall  jets  in 
co— flowing  streams# 
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fi  G.  3-16  Prediction  of  growth  in  the  vertical  direction 


Circular  Wall  Jets  in 
Stagnant  Flow 
(Rajaratnam  &  Pani) 
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FIG.  3-1  7  Prediction  of  growth  in  the  lateral  direction 
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Prediction  of  shear  velocity  at  a  downstream 
centreline  section 


FIG.  3-18 


IV#  Conclusions  and  Recommendations 


In  this  research,  consideration  has  been  made  of  its 
application  to  practical  cases  of  effluent  disposal  in  river 
environments#  The  case  of  one  isolated  diffuser  nozzle  is 
unlikely  to  be  encountered  for  reasons  explained  earlier 
and,  as  previously  shown,  a  common  disposal  technique  is  to 
use  a  series  of  identical  nozzles  with  equal  discharges 
flowing  from  each#  Thus  it  can  be  seen  that  the  situation 
investigated  in  this  work  will  only  exist  for  a  short 
distance  downstream  of  the  nozzle  before  interaction  with 
adjacent  diffuser  flows  occurs# 

The  point  of  occurence  of  this  overlap  may  be  predicted 
from  the  graph  of  transverse  length  scale  growth  for  a 
particular  velocity  ratio  and  the  designer  may  then  choose 
how  far  apart  to  space  the  injection  nozzles  such  that  he 
gets  a  sufficient  reduction  of  concentration  by  the  time 
overlapping  occurs# 

From  then  on  the  jet  could  be  considered  to  act 
approximately  as  if  it  had  been  emitted  from  a  long  narrow 
slot,  although  exact  interaction  of  this  type  has  not  been 
previously  investigated# 

The  results  gathered  from  this  research  and  the 
analysis  applied  to  them  has  shown  that,  as  has  been 
previously  found  for  other  jet  types,  all  sections 
downstream  of  a  jet  have  velocity  profiles  in  the  lateral 
and  free  mixing  regions  that  are  similar#  These  profiles  can 
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be  adequately  described  by  an  exponential 
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relation  given  in 

Chapter  2* 

This  research  has  been  restricted  to  the  measurement  of 
velocity  profiles  in  the  vertical  direction  along  the 
centreline  and  in  one  lateral  plane  at  the  point  of  maximum 
velocity  occurence*  So  the  actual  area  of  investigation, 
compared  to  the  total  size  of  the  Jet,  is  small*  However 
Pani  (1972),  after  studying  4  lateral  sections  at  various 
heights  above  the  bed,  showed  that  all  these  profiles,  no 
matter  whereabouts  in  the  outer  region  they  are  measured, 
can  again  be  described  by  the  exponential  relation*  Thus  any 
point  in  the  jet  can  be  selected  and  its  velocity  predicted 
using  this  relation  and  knowing  the  centreline  velocity  from 
the  vertical  profile  corresponding  to  its  particular  height 
above  the  bed* 

It  has  also  been  shown  that  the  concentration  variation 
downstream  of  a  pollutant  injection  is  of  the  same  form  as 
the  velocity  variation,  both  spread  and  decay  being 
corresponding  functions  of  the  distance  downstream* 

The  justification  for  the  scales  used  in  producing 
non-dimensional  profiles  has  been  shown  theoretically  and 
the  graphs  produced  using  these  scales  have  displayed  gt  S 
agreement  with  previous  research* 

As  was  seen  from  the  graph  of  vertical  jet  growth  a 
levelling  out  occurs  indicating  that  after  a  certain  point 
the  jet  and  main  channel  flow  have  become  one  Inhibiting 
it's  further  expansion*  This  was  also  noticed  on  comparison 
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of  a  series  of  photographs  showing  lateral  spread.  However, 
in  this  case,  it  was  felt  that  the  side  walls  of  the  flume 
were  causing  a  certain  amount  of  constriction  to  growth.  It 
seems  reasonable  that  after  a  certain  time  all  the  excess 
momentum  of  the  jet  will  have  diminished  to  zero  and  from 
then  the  only  diffusion  that  will  take  place  will  be  because 
of  the  ambient  turbulence  in  the  river. 

The  graph  of  lateral  spread  shows  no  indication  of 
le ve lling  off  but  it  only  c  overs  a  range  of  sixty  diameters 
downstream  whereas  the  aforementioned  effects  took  place  at 
least  one  hundred  diameters  away  from  the  nozzle. 

The  actual  amount  of  upward  spread  of  a  pollutant  will 
depend  greatly  on  the  circumstances  of  its  disposal.  In  mo  t 
river  situations  the  effluent  will  be  somewhat  warmer  (and 
therefore  less  dense  )  but  it  will  also  contain  a  certain 
amount  of  dissolved  and  suspended  solids. 

One  area  of  further  research  would  be  an  investigation 
into  the  effects  of  buoyancy  in  the  jet,  due  to  either  heat 
in  the  water  or  perhaos  having  the  jet  inject  into  a  salt 
water  environment.  It  would  be  anticipated  that  the  vertical 
spreading  would  be  significantly  larger,  however  it  is 
doubtful  if  the  lateral  spread  would  be  much  different. 

An  area  that  is  not  well  suited  to  this  type  of  data 
collection  is  the  case  of  a  jet  that  enters  the  surrounding 
channel  flow  with  approximately  the  same  velocity.  Due  to 
the  nature  of  the  apparatus  it  functions  best  when  large 
velocity  variations  have  to  be  detected.  This  means  that  it 
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is  most  suited  to  large  velocity  ratios* 

It  was  noticed  in  Run  no*  3  (  alpha  =  2*21  )  that  t 
effects  diminished  within  such  a  short  distance  from 
potential  core  that  the  data  was  for  some  purposes 
unusuable*  Because  of  the  boundary— layer  shape  of  th 
background  flow  it  was  impossible  to  classify  it  as 
infinite  medium  flowing  past  with  a  constant  average 
velocity*  The  boundary-layer  shape  meant  that  the  ma 
velocity  in  some  cases  occured  near  the  surface,  rat 
in  the  jet* 

When  water  is  Injected  at  a  ratio  of  between  1 
its  effects  will  diminish  In  such  a  short  distance  t 


diffusion  by  channel  turbulence  must  be  considered, 
proposed  to  Investigate  some  spreads  for  a  ratio  of 
1  using  salt  injection  and  a  conductivity  probe  howe 
and  other  considerations  precluded  this* 

The  presence  of  ri ver  cu  rva  ture  would  be  anothe 
important  factor  and  an  investigation  into  the  growt 
jets  in  curved  channels  would  reveal  just  how  much  e 
river  bend  would  have  in  restricting  or  enhancing  di 
of  a  pollutant* 

The  results  gathered  in  this  research  agree  ver 
closely  with  the  predictions  made  by  previous  observ 
the  variation  of  the  jet  as  it  as  it  gets  further  aw 
its  point  of  injection.  The  prediction  of  velocity 
distributions  downstream  of  uniform  circular  wall  je 
as  shown  by  Pani  (  1972  ),  regular  bluff  wall  jets,  in 
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co— flowing  channels  of  large  depth  and  width  to  nozzle  size 
ratio  can  now  be  made  accurately.  From  this  a  suitable  value 
of  pollutant  concentration  can  be  calculated  that  will 
enable  the  effectiveness  of  any  similarity  shaped  effluent 
disposal  system  to  be  determined  without  resorting  to 
expensive  physical  model  studies. 
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Appendix  A. 


List  of  symbols. 


List  of  symbols 


Lenqth  from  position  of  maximum  velocity  to  oosition 
of  half  the  maximum  velocity. 

Channel  width. 

Concentration  of  a  pollutant  at  any  point  downstream 
of  its  injection. 

Maximum  concentration  of  the  pollutant. 

Nozzle  diameter. 

Size  of  bed  roughness  projections. 

Nikuradse's  equivalent  sand  grain  roughness. 

Excess  momentum  of  the  jet. 

Total  Dressure  recorded  by  the  pitot  tube. 

Static  pressure  at  the  wall. 

Nozzle  Reynolds  number. 

Absolute  velocity  at  any  point. 

Absolute  maximum  velocity  at  any  centreline  section. 
Velocity  relative  to  the  channel  at  anv  point. 
Relative  maximum  velocity  at  any  section. 

Initial  velocity  of  the  jet. 

Average  channel  velocity. 

Downstream  or  axial  direction. 

Direction  perpendicular  to  the  wall. 


Transverse  or  lateral  direction. 


✓ 


List  of  symbol s (cont. ) 


Ratio  of  jet  velocity  to  average  channel  velocity. 

Height  from  the  wall  to  the  point  of  maximum 
velocity  occurence. 

Dimensionless  distance  from  the  centreline  of  the  jet. 
Kinematic  viscosity  of  a  fluid. 

Fluid  density. 

Momentum  thickness. 

Shear  stress. 


Shear  stress  at  the  wall. 


Aopendix  B. 


Results  from  Runs  1  to  7. 
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